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THE COSMIC GAMMA-RAY SPECTRUM 
FROM SECONDARY PARTICLE PRODUCTION IN THE METAGALAXY 

F. W. Stecker* 

ABSTRACT 

The purpose of this paper is to discuss the form and intensity of the 
spectrum of cosmic gamma-rays resulting from the production and decay 
of neutral pi-mesons produced in metagalactic cosmic-ray p-p collisions. 
It is assumed that intergalactic space contains ionized hydrogen gas at a 
density of 10“ 5 cm -3 as is consistent with recent X-ray observations at 
0.25 keV. 

Using the Friedmann solution to the Einstein field equations of general 
relativity as adescription of our expanding universe, a discussion is pre- 
sented of the effects of red-shift and spatial curvature on the generation 
and distortion of the local gamma-ray spectrum from the decay of neutral 
pi-mesons. Numerical calculations are presented for the Einstein-de Sitter 
solution, which is found to be an adequate model for these calculations. 
Two models are presented to represent the possible flux of metagalactic 
cosmic-rays. In calculating metagalactic gamma-ray spectra, the effect 
of gamma-ray absorption at large red-shifts is taken into account. 
*NRC-NASA Resident Research Associate 
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A discussion of the results is given. The results indicate that future 
gamma-ray experiments in the 1-100 'MeV region may yield valuable 
information relating to cosmology, cosmogony, and the metagalactic 
cosmic-ray flux. In particular, the metagalactic gamma-ray spectra 
predicted tend to peak near 70 (1 * z ma x)” 1 w ^ ere z max » maximum 
red-shift at which cosmic rays are produced, may correspond to the age 
of the universe at the epoch of galaxy formation. 


THE COSMIC GAMMA-RAY SPECTRUM 
FROM SECONDARY PARTICLE PRODUCTION IN THE METAGALAXY 

INTRODUCTION 


In recent investigations (Stecker, 1967; Stecker, Tsuruta and Fazio, 1968) 
the author made use of recent accelerator and cosmic-ray data to determine the 
details of the cosmic gamma-ray spectrum from the secondary particles pro- 
duced by cosmic-ray collisions in the galaxy. The purpose of this paper is to 
determine the cosmic gamma-ray spectrum from secondary particles produced 
by cosmic-ray collisions in the metagalaxy. This spectrum will differ from the 
galactic (or local) gamma-ray spectrum because moot of the generating collisions 


take place at large distances where we are looking back to a time 


when the uni= 


verse was more compact and collisions were more frequent. These ’’early" 
gamma-rays will be of lower energy due to the progressive red-shift of the gen- 
eral cosmic expansion. Although various estimates of the flux of these meta- 
galactic gamma-rays have been made (Ginzburg and Syrovatskii, 1964a, b; Gould 
and Burbridge, 1965; Garmire and Kraushaar, 1965), none of these workers have 
taken cosmological factors into account in order to properly calculate a spectrum, 
as has been done in calculating Compton X-ray spectra (Cheng, 1967; Brecher and 
Morrison, 1967; Silk, 1968). 
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THE COSMOLOGICAL EQUATIONS 

For the purpose of these calculations, we may consider models of the uni- 
verse which are both homogeneous and isotropic on a large scale. Such models 
can, in general, be described by the Robertson-Walker line element 

ds 2 = c 2 dt 2 - d'C' 2 = c 2 dt - R 2 (t) du 2 . (1) 

Gamma-rays travel along geodesics such that ds 2 = 0. Assuming a gamma-ray 
is emitted at a time t e in an interval At e and received at time t r in an interval 
At r , it can be shown that the relation between red-shift and radius of the universe 
given by 


R(t r ) 

— — — = 1 + z (t ) (2) 

R(t ) ^ e ' K 1 


It follows from (2) that in a universe where most of the energy density is in the 
form of matter 


n (t e ) 

«(*,) 


=(1 + z ) 3 


(3) 




= (1 + z) 


(4) 
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and 


n y(t e ) 

n y (t r ) 


(1 + z) 3 


( 5 ) 


where n (t), T y (t) and n y (t) are the average particle density of matter and 
temperature and photon density of cosmic blackbody radiation in the universe, 
(For a more detailed discussion of the cosmological relations, the reader is 
referred to some excellent articles by Sandage (1961a, b; 1962),) 

We will hereafter designate local (z - 0) quantities with a subscript zero. 

Let G g (E y ) be the gamma-ray spectrum generated by the galactic cosmic-ray 
spectrum, I g (E p ) in traveling a unit particle length (1 cm" 1 ) through the inter- 
galactic medium. (This spectrum is the same as the quantity I (E y )/<nL)> cal- 
culated by Stecker (1967) and Stecker, et. al. (1968).) 

We now assume that some ubiquitous generating mechanism causes cosmic- 
rays to be produced with the same power law throughout the universe as observed 
at the earth, so that the metagalactic cosmic-ray spectrum differs only .in ab- 
solute intensify from the galactic cosmic-ray spectrum. It follows that the form 
of the cosmic gamma-ray spectrum anywhere in the metagalaxy, when observed 
in the co- moving frame at that point, will be the same as the form of G g (E r ). 

We may then write down an expression for the integrated metagalactic gamma- 
ray flux in any direction as 
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( 6 ) 


I (E y) ^ 



din (tolfl 


Q,(M> e -r(E 7 ,4) 

(1 + z (4)) 4 e 


where the factor, (1 +z ) 4 , takes into account the reduction in flux due to the time 
dilation and volume diminuation factor and e“ r represents absorption of gamma- 
rays along the line of sight, I g is the galactic cosmic-ray flux and 1{1) is the 
cosmic-ray flux at a distance l . Equation (6) may be put into a much more 
convenient form by expressing it as an integral over z . We then obtain 




dz n (z) 


I(z) 


f(E y , z) ~r (E y , z) d'C 
(1 + z) 4 dz 


( 7 ) 


Since the energy of a gamma-ray is directly proportional to its frequency, it 
follows that 


G g (E r ,z) = (1 + z) G g [(1 + z) E y ] 
It also follows from (3) that 

n (z) =n 0 (l + z) 3 


The quantity 


dl 

dz 



( 8 ) 


(9) 


(10) 
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depends, in general, both upon :.ne cosmological model involved and the epoch of 
world-time which defines the acceleration (or deceleration) of the expansion. 

In Friedmann-type solutions to the Einstein equations, it is found that the ex- 
pansion of the universe is decelerating. The magnitude of this deceleration is 
usually denoted by the deceleration parameter q . In the usual notation, the Hubble 
expansion parameter, H, and the quantity q are defined by the relations 


H 


q £g) 

R(t) 


( 11 ) 


and 


q Q 


g (t) 

R (t) H ;? 


(12) 


In a decelerating universe, therefore, q > 0. 

For the Einstein-de Sitter model, R (t ) can be expressed explicitly in terms 
of t by the relation 


R (t) = (677G/PR 3 ) 173 t 2/3 


From (11) and (13) it then follows that 

q = i f or all t , 
2 


(13) 


(14) 
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The Einstein-de Sitter model is a good approximation to the universe if it 
has not yet reached a highly evolved state. It is also eompatable with the most 
probable values of q as discussed by Sandage (1961a, 1962), based on the observed 
magnitude -red shift relation, and is consistent with the X-ray observation by 
Henry, Fritz, Meekins, Friedman and Byram (1968) interpreted to be brems- 
strahlung from a metagalactic gas having a density af the order of 10“ 5 cm" 3 , 

Under the assumption of a Euclidean model, we will now determine the 
cosmological effects on the metagalactic gamma-ray spectrum. It can be shown 
(Sandage, 1961b) that 


d* . cHq 1 

dz (1 + z ) 2 (1 + 2 q 0 z ) 1 ' 2 


(15) 


where cH" 1 = 10 28 cm. In the Einstein-de Sitter case, q 0 = 1/2 and we may 
take in equation (7) 


= 10 28 .. ( 16 ) 
dz (1 + 2 ) 5/2 

ABSORPTION OF METAGALACTIC GAMMA-RAYS 

An excellent discussion of the absorption processes affecting cosmic gamma- 
rays has been given by Fazio (1967). The principal absorption process to be 
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considered is that of electron-positron pair production through interaction with 
the universal black-body radiation field, i.e., the reaction 

y + y -* e + + e” (17) 

Detailed calculations of the energy-dependent absorption probability for 
this process have been performed by Gould and Schreder (1967). They have 
shown that for a gamma-ray of energy interacting with a black-body radiation 
field of temperature T y 


dr ^ a 2 
dl ” 2'jtK 



SI e-* 


(18) 


where 


p = ( m ° 2 ) 2 » i 

kI y E y 

a - 1/137 is the fine -structure constant, A -h /me - 3.86 x 10" 11 cm, and k 
here is Boltzmamds constant. The local black-body temperature has been found 
by Stokes, Partridge and Wilkinson (1967) to be 

T 0 = 2.7 °K. (19) 

so that the condition g » 1 corresponds to the condition 
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1.12 x 10 6 GeV 


( 20 ) 


R y « 


(11 + z) 2 


(see equation (4) and (8)). 

We will restrict ourselves here to a determination of the gamma-ray spectrum 
below 1 GeV and z £ 10 2 (as will be discussed later) so that the approximation 
given by equation (18) will be generally valid. Therefore, from (16) and (18), 
we find* 


r(E y ,z) = 3.9 x 10 8 E- 1 / 2 



dy 


~1.12 x 10 6 ~ 

q + y) 2 E y ] 
(1 + y) 1/2 


(21) 


THE METAGALACTIC COSMIC-RAY SPECTRUM 

It now remains only to specify a suitable model for the metagalactic cosmic 
ray flux. We will assume that at some early epoch, corresponding to z £ z max 
conditions were unsuitable for the acceleration of cosmic-rays. We will con- 
sider z max to correspond to the epoch of galaxy formation and consider two 
possible models for the origin of a metagalactic cosmic-ray flux. For model I , 


*As an intermediate solution to the problem considered in the text, numerical solutions were ob- 
tained for the implicit relation, r( E y , z y ) = 1, which defines the red-shift, z y , beyopd which the 
universe becomes opaque to gamma-rays of local energy E y . It was found that the numerical solu- 
tion to equation (21) may be quite well approximated by the relation 1 + z - 2.30 x 10 2 E y °' 484 
with E y expressed in GeV. Thus*, although absorption was taken into account in this calculation, 
the effect is negligible as the universe may be considered to be transparent out to z y . 
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we will assume that the metagalaotic flux arises through leakage from the halos 
of radio galaxies from z = z max to z - 0. For model IX, we will assume that 
this flux was created primarily in a burst at the time of galaxy formation. (Other 
possibilities will be considered in a future paper.) For z max , we will also con- 
sider two extremes. One extreme is z max = 10 3 , which corresponds to the 
earliest epoch when galaxy formation could probably occur. At z = 10 3 , the 
black-body temperature of the universe was of the order of 10 3 - 10 4o K, cool 
enough for ionized hydrogen to combine to form a neutral gas. According to 
Peebles (1965), z ~ 10 3 also corresponds to the epoch when gas clouds may 
begin to form gravitationally bound systems. 

The other extreme for z max which we may consider corresponds to the 
highest red“Shift yet observed for a quasar, vis., 2.2. This is, of course, an 
extreme which is technique- limited rather than being limited by any physical 
criteria, and it is included mainly for purposes of discussion. We will also 
consider various intermediate values for z max of 4, 9, and 10 2 . (Doroshkevich, 
et. al. (1967) suggest that galaxy formation took place at z = 10 - 20 whereas 
Weymann (1967) suggests z = 10 2 .) 

It is important to note here that the upper limit, z max , may be effectively 
restricted, not by the epoch of galaxy formation, but by attenuation of the meta- 
galactic cosmic-ray flux due to the collisions themselves. The cross-section 
for inelastic cosmic-ray p-p collisions is of the order of 30 mb. Therefore, 
the lifetime of the metagalactic cosmic-rays against collisional losses is given by 
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r c = 


ncrc 


n Q crc (1 + z)‘ 


- 10 ?0 (1 + z)" 3 sec 


( 22 ) 


for n 0 = 10" 5 cm" 3 . 


The lifetime of the universe at a red-shift z is given by 


7* u - 10 17 (1 + z)" 3//3 sec 


(23) 


Cosmic-rays cannot accumulate in the metagalaxy if the ratio, t c /t u < 1. 
The condition r c /r u = 1 therefore defines a critical value of z max = 10 2 beyond 
which a further buildup of metagalactic cosmic-rays cannot occur. With these 
limitations on z in mind, we will now consider the various ideal models for 
describing the metagalactic cosmic-ray flux. 

For model I, we assume a constant leakage rate so that the total number of 
cosmic-rays in the metagalaxy is proportional to the time elapsed since galaxy 
formation. It follows from (13) that this time is given by 


r g S 10 10 [(1 + z)- 3/2 - (1 + z n , ax )- 3/2 ] yrs. (24) 


The cosmic-ray density will then increase with red-shift according to the 
relation 
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~ (1 + z) 3 C(1 + z)" 3/2 - (1 + z max )‘" 3/2 3 


(25) 


X'(g) 


However, the cosmic-ray^ which produce the neutral pi-mesons necessary for 
gamma-ray production are only those above a threshold kinetic energy, E th -M , 
of about 300 MeV (Stecker, 1966). We must therefore determine 


I 00 = I # (E>E th ; z) 


(26) 


For a power law cosmic-ray spectrum of the form 

I(>E)~E- 3/2 , (27) 

as is approximately valid in the energy region where most pi-mesons are 
produced (Stecker 1967), it follows from the red-shift relation that 

I (&> E th ; .) = !' U>~^\ 


= r (z) 


1 + z 


1 + Z 


max 


3/2 


(28) 


so that we must use an effective flux of 


— ~ (1 + z) 3 


1 + z 


1 + z 


max J 


3/2 


C(1 + z)“ 3/2 - (1 + z max )" 3/2 ] 
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For model II, we assume that the metagalactie cosmic-rays were created 
in a burst at the time of galaxy formation. We thus find for this model that 


~ (1 + z) 3 

B 


1 + z 


*1 3/2 


1 + Z 


max 


(30) 


Using the models defined by equations (29) and (30) , together with equations 
(7), (8), (9), (16), and (21), the metagalactie gamma-ray spectra produced by 
models I and II were calculated. The results were normalized by requiring the 
integral gamma-ray spectrum above 100 MeV to be equal to 1.1 x 10" 4 cm -2 sec"* 1 
sr _1 , according to the results of Clark, Garmire and Kraushaar (1968), measured 
by the detector aboard the OSO-3 satellite. A 3 has been noted previously, such 
a normalization makes possible the determination of upper limits on the value 
I 0 /I g , i.e., the present metagalactie intensity of cosmic-ray nucleons. The upper 
limits are given by Stecker (1968). 

The gamma-ray fluxes thus calculated are given in Figures 1-4. Figure 5 
shows the gamma ray flux expected from the. galactic halo in the direction of the 
pole, taking <(nL)> = 3 x 10 20 cm“ 2 and based on previous calculations (Stecker, 
1967; Stecker, et. al. 1968).* It can be seen that the local gamma-ray flux from 
the galactic halo will not explain the data of Clark, et. al. and should be unimportant 
compared to the metagalactie flux. The metagalactie gamma-ray spectra tend to 
* Thus Ipol. (E r ) = 3xl0 2 °G g (E r ). 
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to peak near 7 x 10" 2 /(1 + z max ), GeV, being weighted toward higher red-shifts 
by the effect of greater densities at earlier epochs. Because of the density effect, 
a cosmic-ray burst at large red-shifts is much more effective in producing gamma- 
rays then a continuous production of the same number of cosmic-rays. 

CONCLUSIONS 

Present evidence about the flux of cosmic-rays between the galaxies is quite 
meager. The most promising way to study the flux is by a satellite experiment 
measuring the isotropic gamma-ray flux in the region between 1 and 100 MeV. 

Such gamma-rays can supply us with direct information on metagalactic cosmic- 
rays, because they travel to us in straight lines and suffer little absorption. 
Theoretical metagalactic gamma-ray fluxes from tt° decay are presented here 
under various assumptions as to the metagalactic cosmic-ray flux. These pre- 
dictions indicate that an experimental determination of the isotropic gamma-ray 
spectrum at high galactic latitudes and in the energy range 1 - 100 MeV, could 
supply valuable information, not only about metagalactic cosmic-rays, but also 
about such fundamental questions as when the galaxies were formed, since the 
metagalactic gamma-ray spectrum will peak near 70 (1 + z max ) _1 MeV. 

ACKNOWLEDGMENTS 

The author is indebted to Dr. Richard C. Henry of the Naval Research Lab- 
oratory for supplying him with the details of the NRL X-ray observations prior to 
publication and for several stimulating discussions. The author would also like 


13 



to acknowledge and thank Dr. Frank C. Jones and Dr. Joseph Silk for their com- 
ments and Mr. Joseph Bredekamp for programming the numerical calculations 
essential to this paper. 


14 


REFERENCES 


Brecher, K. and Morrison, P., 1967, Astrophys. Jour. Letters 153, L61. 

Cheng, C. C., 1967, Nature 215, 1035. 

Clark, G. W., Garmire, G. P. and Kraushaar, W. L., 1968, Astrophys. Jour. 
Letters 153 , L203. 

Doroshkevich, A. G., Zel’dovich, Y. B., and Novikov, I. D., 1967, Sov. Ast. A. J., 
11, 233. 

Fazio, G. G., 1967, Ann. Rev. Astron. and Ap., 5, 481. 

Garmire, G. and Kraushaar, W., 1965, Space Sci. Rev., 4, 123. 

Ginzburg, V. L. and Syrovatskii, S. 1., 1964a, Sov. Phys. JEPT, 18, 245. 

Ginzburg, V. L. and Syrovatskii, S. X., 1964b, The Origin of Cosmic Rays , Mac- 
millan Co., New York. 

Gould, R. J. and Burbidge, G. R., 1965, Ann. d’ Astrophys., 28, 171, 

Gould, R. j. and Schreder, G. P., 1967, Phys. Rev., 155, 1404. 

Henry, R. C., Fritz, G., Meekins, J. F., Friedman, H., and Byram, E. T., 1968, 
Ap. J. Letters 153 , Lll. 

Peebles, P. J. E., 1965, 142, 1317. 


15 


Sandage, A., 1961a, Ap. J., 133, 355, 
bandage, A., 1961b, Ap. J., 134 , 916. 

Sandage, A., 1962, J. Soc. Indus, Appl. Math, 10, 781. 

'Silk, J,, 1968, Ap. J. Letters 151, L19. 

Stecker, F. W., 1966, Smithsonian Ap. Ohs. Spec. Kept. No. 220. 

Stecker, F. W., 1967, Smithsonian Ap. Obs. Spec. Kept. No. 260. 

Stecker, F. W., Tsuruta, S., and Fazio, G. G., 1968, Ap. J., 151, 881. 

Stecker, F. W., 1968, Nature, In the press. 

Stokes, G. A., Partridge, E. B., and Wilkinson, D. T., 1967, Phys. Rev. Letters, 
19, 1199. 

Weymann, R., 1967, Ap. J., 147, 887. 


16 



FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


. Metagalaotie differential gamma-ray spectra from cosmic-ray p-p 
interactions based on a cosmic-ray flux produced by constant leakage 
from radio galaxies (Model I) and shown for various maximum red- 
shifts as discussed in the text. 

. Metagalactic integral gamma-ray spectra from cosmic-ray p-p inter- 
actions based on a cosmic-ray flux produced by constant leakage fe orn 
radio galaxies (Model I) and shown for various maximum red-shifts 
as discussed in the text. 

. Metagalactic differential gamma-ray spectra from cosmic-ray p-p 
interactions based on a cosmic-ray flux produced by a burst of cosmic 
rays at z max (Model II) as discussed in the text. 

. Metagalactic integral gamma-ray spectra from cosmic-ray p-p inter- 
actions based on a cosmic-ray flux produced by a burst of cosmic 
rays at z mQX (Model II) as discussed in the text. 

. The expected gamma-ray flux from the galactic halo in the direction 
of the pole, taking <(nL)> = 3 x 10 20 cm" 2 and based on previous calcu- 
lations (Stecker, 1967; Stecker, et. al, 1968). 
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Figure 2-lntegral Spectrum for Leakage Model (Model I) 
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